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Abstract 

We study the shape differentiabihty of a cost function for the flow of an incom- 
pressible viscous fluid of power-law type. The fluid is confined to a bounded 
planar domain surrounding an obstacle. For smooth perturbations of the shape 
of the obstacle we express the shape gradient of the cost function which can be 
subsequently used to improve the initial design. 

Keywords: shape optimization, shape gradient, incompressible viscous fluid, 
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1. Introduction 

Shape sensitivity analysis is a fundamental step towards the numerical solu- 
tion of the shape optimization problems governed by partial differential equa- 
tions. In the context of Navier-Stokes equations, the sensitivity analysis of shape 
functionals is performed in [lil Q for the incompressible case, and in [l7j for 
the compressible case. Usually, in numerical solution of applied problems this 
step is made formally, see e.g., for the related results. 

The present paper is focused on the rigorous analysis and new results for a 
nonlinear nonstationary model in two spatial dimensions are proved. We refer 
the reader to ^ for some results in the case of a stationary model. 

We consider the time-dependent flow of an incompressible fluid in a bounded 
domain fl -.^ B \ S in K^, where B is a container and S is an obstacle. Motion 
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of the fluid is described by the system of equations 



) + Cv 


= f 


in Q, 


divv 


= 
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V 


= 


on E, 


v(0,.) 


= Vo 


in ri. 



[PiSi)) 



Here Q := (0,T) x fi, E := (0,T) x 917, where (0,T) is a time interval of 
arbitrary length, v, p, C, f stands for the velocity, the pressure, the constant 
skew-symmetric Coriolis tensor and the body force, respectively. The traceless 
part § of the Cauchy stress can depend on the symmetric part Dv of the velocity 
gradient in the following way: 

§(Dv) = i/(|Dv|2)I[])v, (1) 

where V, IBvp is the viscosity and the shear rate, respectively. In particular, we 
assume that v has a polynomial growth (see Section [53] below), which includes 
e.g. the Carreau and the power-law model. 

In the model the term of Coriolis type is present. This term appears e.g. 
when the change of variables is performed in order to take into account the flight 
scenario of the obstacle in the fluid or gas. 

The aim of this paper is to investigate differentiability of a shape functional 
depending on the solution to ( |P(11)| with respect to the variations of the shape 
of the obstacle. We consider a model problem with the drag functional 



J(n) := / / (S(Dv) - pT)n ■ d, (2) 
Jo Jas 

with a given constant unit vector d. Instead of J one could take other type of 
functional, since our method does not rely on its specific form. 

Fluids whose viscosity depends on the shear rate through ^ form an impor- 
tant class of non-Newtonian fluids (see e.g 21, 2^ [3] for general references). 



Particular models find applications in many areas of chemistry, biology, medicine 
or engineering. The first mathematical results were established already 40 years 
ago in [n, 12, 13 1, for recent references see e.g. 0, H, IH, Sj. In context of op- 



timal control the fluids with shear dependent viscosity were studied in |22l. |27|. 
Numerical results for a shape optimization problem can be found in [ij , see also 

Our main interest is the rigorous analysis of the shape differentiability for 
( |P(1I)[ ) and ([2]) . We follow the general framework developed by Sokolowski and 
Zolesio using the speed method and the notion of the material derivative. 
Let us point out that due to ([T]) the state problem is nonlinear in its nature. 
In such cases the shape sensitivity analysis is usually restricted to 'small data', 
i.e. for sufficiently small body forces, initial conditions or short times. For the 
presented results no such restriction is necessary, because of uniqueness and 
regularity of the state variables. 
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1.1. Shape derivatives 

We start by the description of the framework for the shape sensitivity anal- 
ysis. For this reason, we introduce a vector field T g C^(E^,E^) vanishing in 
the vicinity of dB and define the mapping 

y(x) =x + eT(x). 

For small e > the mapping x h-> y(x) takes diffeomorphically the region Q, onto 
n^=B\Se where = y{S). Denoting :^ (0, T) x n^, := (0, T) x dn,, 
we consider the counterpart of problem ( |P(r^)[ ) in Q^, with the data f|Q^ and 
Vqio^- The new problem will be denoted by {P{il^)) and its solution by (v£,pe). 

For the nonlinear system (P{Q)) we introduce the shape derivatives of solu- 
tions. To this end we need the linearized system of the form: 

Find the couple (u, tt) such that 



dtu + div [u (g) V + V (g) u - S'(Dv)]ri)u] + 



Vtt + Cu 


= F 


in Q, 


div u 


= 


in Q, 


u 


= h 


on S, 


u(0,-) 




in f2, 



(Plinm) 



where F, h and Uq are given elements. 

The shape derivative v' and the material derivative v of solutions are for- 
mally introduced by 

V := Imi , V lim . 

e^O £ e-i-0 £ 

The standard calculus for differentiating with respect to shape yields that v' is 
the solution of (Piin(i^) ) with the data F = 0, Uq = 0, and h = — (?v/9n(T • n). 
Using p3p as the definition of J we obtain the expression for the shape gradient: 

dJ(f7;T):=lim:^(M-M 

e^O £ 

^ I ^'{T)-i + [ [(Cv')-^ + (§'(rov)Dv'-v'«)v-v®v') : V|] 
Jn Jq 

- r f (f.d)T.n. (3) 

Jo Jas 

In the above formula, the part containing v' depends implicitly on the direction 
T. This is not convenient for practical use, hence we introduce the adjoint 
problem for further simplification of (|3|): 
Find the couple (w, s) such that 



-9tw - 2{Bw)v - div [§'(Dv)^Dw] + Vs - Cw = 

div w = 



in Q, 

w = d onS, (Padj(f^)) 
w(T, •) = m n. 
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Consequently, the expression for dJ reduces to 



dJ{n;T) 



OS 



(§'(Dv)TDw - sI) : ^ ® n + f • d 
^ ' an 



T n. 



(4) 



In order to prove the result given by ^ and Q we need the material deriva- 
tives. In particular, it is sufficient to show that the linear mapping 

T ^ (iJ(0;T) 

is continuous in an appropriate topology, see the structure Theorem in the book 
23[ for details. 

The paper is organized as follows: In Section [2] we impose the structural 
assumptions and collect the basic facts about the existence of weak solutions to 
^P(iTj\ . The main results of the paper are stated in Section [3l Well-posedness 
of the linearized systems is studied in SectionlH The rest is devoted to the proof 
of the main results. In Section [5] the problem is formulated in the fixed domain. 
Section |6] deals with the shape stability. Section [7] with the Lipschitz estimates 
and the existence of the material derivative and Section |5] with the existence of 
the shape gradient of J. 



2. Preliminaries 

We impose the structural assumptions on the data, state the known results 
on well-posedness of ( |P(r2)| and introduce the elementary notation for shape 
sensitivity analysis. 

2.1. Structural assumptions 

We require that S has a potential $ : [0, oo) — > [0, oo), i.e. §y (B) = 
93>(|Dp)/aDi-,-. Further we assume that $ is a function with $(0) = 
and that there exist constants Ci , C2 , C3 > and r > 2 such that 

Ci(l + |Ar-2)|B|2 <§'(A) :: (B ® B) < C2(l + |A^-2)|B|^ (5a) 

|S"(A)| <C3(l + |Ar-3) (5b) 

for any 7^ A,B e ^lym- The above inequalities imply the following properties 
of S: 

Lemma 1. LetS satisfy ([5]). 

(i) There is a constant C4 > such that for every A G Rsym-' 

|§(A)| <C4(l + |Ar-i). (6) 

(ii) There is a constant C5 > such that 

(§(A) - §(B)) : (A - B) > C5 |A - Bf (7) 
for all A, B G R^x^. 
The proof can be found e.g. in [l5l |. 
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2.2. Weak formulation 

For the definition of tlie weak solution we will use the space 



Wj;L(f^) :={0e Wj''-(f}); dive/.: 



Let vo e Wo;;ji^(rj) and f e L^(0, T; (Wg;^; jrj))*). Then a function v € 

L'-(0,T;W^;;;iJ17)) with ^tv e ^^(o, T; (W^;^;^ (17))*) is said to be a weak 
solution to the problem (P(0)) if v(0) = vq and 

nT 



0}. 



§(Dv) : - V ® V : V(/> + Cv • 



= / f • (8) 
Jo 



for every e L'^(0,T;W;';;i^(n)). Note that the pressure is eliminated since 
test functions are divergence free. 

The existence and uniqueness of a weak solution to (P(il)) with C = in two 
space dimensions is due to Ladyzhcnskaya [H, 12| and Lions [l^l- Regularity 
of weak solutions was studied e.g. by Malek et al. [3], Beirao da Veiga et al. 
ii, Kaplicky 0. We recall the following result: 



Theorem 2 (Kaplicky [l3|)- ^et C &e a bounded domain with C^+^ 
boundary and T > 0. Let ([5]) hold for some r £ [2,4), C = and 

f eL°°(0,T;L2(r!)), i{0)eL^{n), dtieL'{0,T;{Wl%jn)r), (9) 

voeWP''{n)nwl%^{n), p > 2 (p = 2 */r = 2). (10) 

Then the unique weak solution v, p of {P{il)) satisfies for s G (1,2) (s ^ 2 if 
r = 2): 

V e L°°(0, T; W2.-(l])), Stv e L°°(0, T; L2(r!)) n £^(0, T; W^J;^;^ (17)). 
Moreover, if there is a q > 2 such that 

f e L°°(0,T;L'?(17)), e L«(0,r;W-i'«(17)), (11) 

i/ien t/iere exists q > 2 and a > such that for all e £ (0,r) it holds: 
Vv,pe L°°(e,r;W^i'«(l]))nC°^"([e,r] xH). 



We wiU also need the L'' theory for the generalized Stokes system (Piin(f2) ). 
Here and in what follows we will assume that the non-homogeneous boundary 
condition is time-independent, i.e. h = h(x). Theorem 4.1 in Q implies (see 
also (H): 
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Theorem 3 (Bothe and Priiss 0). Let Q e C^'i, T > 0, g G (1, oo) \ {§,3}, 
Dv e C(0), F e L«(Q), h e W2-i/9^9(ar2), and uo G W2-2/9,9(f7)_ i^t h anrf 
Uq satisfy the compatibility conditions: div Uq = 0, h • n = Uq ■ n on dn and for 
g > I h = Uo on dfl. 



Then (Pn^ijl)) has a unique solution (u, tt) in the class 



u e W^i^'?(0, T; L'?(17)) n £"(0, T; W^''^{n)), tt G i«(0, T; W^''^{n)). 

The additional Coriolis term Cv presents a minor technical obstacle in the 
existence analysis: Since C S ^Ikew^ ^ priori estimate remains unchanged, 
namely every weak solution v satisfies the energy equality 

^l|v(t)||^,o+ f I S(ID)v):rov= f f f -v + illvoll^^o, for a.a. t G (0, T). 
^ Jo Jn Jo Jn ^ 

(12) 

Considering a sequence {v^} of approximate solutions that satisfy ([T2|). the 
limit passage 

Cv^ ^ Cv in L'^(0,T; Wi^''(17)) 

is a straightforward consequence of strong monotonicity of § and the weak com- 
pactness of {v^} in L''(0, T; W^'''(r2)). In other aspects the proof of the exis- 
tence and the uniqueness of a weak solution follows line by line the original one. 
Concerning regularity, if we put Cv to the right hand side, the assumptions of 
Theorem [3] are still satisfied, hence the same result holds. 

In order to identify the material derivative of v and to apply Theorem [3] 
we will need Vv G C{Q). Theorem [51 however, guarantees this only for r = 2 
or locally in time. To overcome this technical obstacle, we use the idea of 
Wachsmuth and Roubi'cek [i^] and impose the following restriction on the initial 
condition and the body force. 

Assumption 1. In what follows we suppose that is a bounded domain in 
with C^'^ boundary and r G [2,4), the function f is extended by in the 
complement of Q, and suppvg H dS = 0. Further we assume: 
For r = 2: 

The initial value Vq and the body force f satisfy ([9]), (|10[) and (jlip. 
For r > 2: 

There exists t > Q andi <E i°°(-T, T; L«(17)) withdfi G L«(-t, T; W-i'«(17)) 
such that f = f Iq; the initial value vo is equal to v(0) with v being the solution 
of the problem 

dtv + div (v «) v) - div S(]D)v) + Vp + Cv = f , 

divv = in (— r, 0) x fi, 

v = on{-T,0)xdn, 

v(— T, •) ~ v_T- in il, 

where v_t- satisfies (|10|) with v_t- in place o/vq. 
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Corollary 4. Let Assumption{^hold true. Then there is a unique weak solution 
V to ( |P(r^)[ ) and the associated pressure p which satisfies: 

V e i°°(0,T;W2-''(17)), 

VveC°'"(Q), 

pe L°°(0,T;W^i^9(17))nC(Q) 
/or soTTie q > 2 and a > 0. 

Remark. The above result applies only to the unperturbed domain, i.e. e = 0. 
For the other cases we assume that the solution (v£,pe) belongs just to the class 
of the weak solutions. 

Let us point out that equation ^ which defines J is not suitable for weak 
solutions in general, since the energy inequality does not provide enough in- 
formation about the trace of p and Bv. We therefore introduce an alternative 
definition that requires less regularity. Let us fix an arbitrary divergence free 
function ^ e C'^{B,M?) such that ^ = d in a vicinity of S. Then, integrating 
([2]) by parts and using (P[VL) \ yields: 

J{^1) = / (v(r) - vo) • I + / [(Cv - f) • ^ + (§(Dv) - V ® v) : V|] . (13) 
Jn JQ 

Note that this identity is finite for any v e L'^{0, T; W^-'^{VL)) and v(r) e l7{n). 

2.3. Deformation of the shape 

Arguing similarly as in Section [2?2l we find that (P(J7e)) has a unique weak 
solution V, G U{Q,T-Wl^;^^{n,)) with dt^, € (l^Q^T-W^^l.^in,)))* which 
satisfies the energy inequality 

1„_ , /■* f ^_ /■*/■„_ 1 



:^l|ve(<)ll2.o. + / / §(Bv,):Bv,= / / f •v, + -l|volllo. (14) 

for a.a. t e (0, T) and for t = T. 

Let us introduce the following notation: We will denote by D T the Jacobian 
matrix whose components are (DT)jj = (VT)ji = diTj. Further, 

N(x) := 0(x)M~i(x), M(x) := I + eDT(x), 0(x) := det M(x). 

One can easily check that the matrix N and the determinant q admit the ex- 
pansions: 

fl = l + edivT + 0(e2), N = I + eN' + 0(£^), N' = (div T)I - D T, (15) 
where the symbol 0(£^) denotes a function whose norm in C^{fl) is bounded by 
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The value of the shape functional for is given by 
J{n,) / (v,(T) - vo)-^,+ f [(Cv, - f) • + iH^^e) - V, ® V,) : V| J , 

where := (N^^^) o yj^. Using the properties of the Piola transform (see 
e.g. Theorem 1.7-1 in \^) one can check that div^^ = 0. If and were 
sufficiently smooth, it would hold that 

J{n,)= [ [ {E>{m,) - pel)n ■ d. (16) 

Nevertheless, as opposed to ( |P(r2)| , we do not require any additional regularity 
of the solution to the perturbed problem {P{ils)) and hence the expression in 
(fT6)) need not be well defined. 

We introduce the auxiliary function v: 

N^Ve o y — V 
V := lim , 

e^O £ 

which is related to the material derivative v by the identity 

V = N'^v + V. 

For the justification of the results of the paper we will use v since, unlike the 
material derivative, it preserves the divergence free condition. 

3. Main results 

The first step is the existence of the function v and hence also of the material 
derivative. 

Theorem 5. Let Assumptions^ be satisfied. Then the function v exists and is 



the unique weak solution of (Piin(r2)) with the data 



F = A[, := (N' + N'^ - I tr N')9tv + div (v ® N'^v) + N'div (v ® v) 
+ div [§'(]D)v) (((N' - ItrN')Vv),.j,,„ - ©(N'^v)) + N'^§(]D>v)] 
- N'div §(Dv) + ((N' - I tr N')C + CN'^) v + (I tr N' - N')f + (Vf )T, (17a) 

h = 0, (17b) 
uo = vo N'^vo + (Vvo)T. (17c) 

The following estimate holds: 

l|v||L-=(0,T;L2(n))nL2(0,T;Wi.2(n)) < C'llAo|1^2(o.T;W^'^. (H))' ^ C*]] T|l ^2 (H) ■ 

(18) 
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Due to the term (Vf)T the right hand side Ag is not integrable, thus one 
cannot apply Theorem [3] Instead, the well-posedness will be investigated in 
Sectional The next result concerns the existence of the shape gradient. 

Theorem 6. Under the assumptions of Theorem \5l the shape gradient of J 
reads 

dj{n,T) = Mx) + j,{T), 

where the dynamical part and the geometrical part is given by 

J^(v)= [ (v(r)-vo)-|+ /" [(Cv) •! + (§' (ID)v)Dv-v®v-v (g)v) : V|], 
Jn JQ 



Je(T)= / (ltrN'-N'-N'^)(v(T)-vo)-| 
Jn 

+ J I [ (I tr N' - N') Cv - CN'^ V - (I tr N' - N') f - (Vf )T] • | 

+ N'^v + §'(©v) ((N'Vv - W{N'^^r))syrn - (trN')]D)v) + N'^§(Dv)] ; V| 

+ [v (g) V - S(]D)v)] : V(N'^,e)}, 
respectively. In particular, as v depends continuously on T, the mapping 

dJ{n,T) 

is a bounded linear functional on C^(R^,R^). 

Based on the previous result we can deduce that the shape gradient has the 
form of a distribution supported on the boundary of the obstacle. Since this 
representation is unique, the formal results derived in Section 11.11 are justified 
provided that the shape derivatives and adjoints exist and are sufficiently regu- 
lar. This issue will be addressed in the following section. At this point we state 
the final result. 

Corollary 7. Let Assumption{l\ be satisfied. Then 

(i) the shape derivative v' exists and is the unique weak solution to (Piin(51) I 

with F = 0. h = -T^(T • n), uq = 0; 
on 

(ii) the shape gradient of J satisfies 



(Hi) the adjoint problem (Padj(f^)l has a unique weak solution that satisfies for 
arbitrary 5 G (0,T); 

w e ^2(0, T - (5; W2'2(S1)), s G 2.2(0, T - (5; W'^'^{n)). 

(iv) Finally, if f G £^(0, T; W^'2(r2)) then the shape gradient of J has the 
representation (j4]) in the following sense: 



dJ(n; T) = - lim /" / 
-5^0 Jo Jds 



(§'(Dv)^Dw - si) : — ® n + f • d 
^ ' on 



T n. 

(19) 
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Let us note that the formula ^ is here replaced by ([TO]). This is due to 
the fact that the terminal and boundary conditions for w do not satisfy the 
hypothesis of Theorem |31 In fact, (j?]) holds if we replace the values Dw(T, •) 
and s{T, •) by their L^(9f2)-limits for t ^ T, see LemmalHl 



4. Well-posedness of and of (Padj(f^)l 



In this section wo show that the linearized system has a unique 

solution under assumptions that are weaker than those of Theorem |3l Indeed, 
we have the following result. 

rl.2 



Lemma 8. Let Q G C"'S v e W^'°°{Q), F e 2.2(0, T; W^;^;^ (17)*), h € 

L2(0,TVW^r!)), dth e LHO,T;Wl%jnr) and Uq - h(0, •') S Lg ( 
Then (Pii^{fl)) admits a unique weak solution which satisfies the estimate: 



ll'^llL°'(o,T;L2(n))nL2(o,r;W^'^j^ (O)) - C'(l|uol|2,n + I|F1|^2(o^t;W^^J'^j^ (a)*) 



|h||L2(0,T;Wi-2(J2)) + ||9th| 



Proof. Using the standard Galerkin method, one checks that the approximate 
solutions satisfy 

l|u^llL°°(0,T;L2(n))nL2(0,T;Wi.2(n)) < C, 

where the constant on the r.h.s. depends on ||v|ji^(X)_Q and the respective norms 
of Uq, h, F and |C|. Hence there is a weak limit u of {u^} which is a weak 
solution to (Piin(r2)|. The uniqueness can be proved testing by the difference 
ui — U2 of two solutions, from which one gets for a. a. t € (0,r): 

\\uiit) - U2(<)!||n + - U2)||^,Q < C||v||oo,q||ui - U2||2,q||V(ui - U2)||2,Q. 

Gronwall's identity then directly implies ui = U2. □ 

A direct consequence of the lemma is the existence and uniqueness of the 
shape derivative v'. 

Applying the same technique as in the previous lemma and additionally 
Theorem [3l one obtains the result for the adjoint problem: 

Lemma 9. Let fl G C*^'^, ^ G C°°{B). Then there is a unique weak solution 
{w, s) to ( |Padj(^)| . // in addition il e C^-^ , then for any S e (0, T); 

w e L2(0,T-(5;W2'2(f7)), dtw e L^{0,T-6;L'^{n)), s e L^{0,T~6;W^''^{n)). 

In particular, \ims\o^w{T — S, •) and lim5\,o s{T — S, •) in L^{dVl) exist and 
are finite. 
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Proof. Using the change of variables t i-^ T — t we transform (Padj(i^)l to 



(Fiin(f^) I and obtain uniqueness of the adjoint state w. It is however not pos- 
sible to apply Theorem [3] since the terminal condition for w is nonzero on dS. 
Nevertheless, considering the adjoint problem on the time interval (0, T — S), we 
find that the terminal condition w(r — J, •) is compatible and hence 

!|ID'w!|2,(o,T-5)xan + \\s\\2.{o,T-5)xan < C\d\. 

Existence of the limits of Dw and s follows from the continuous dependence on 
the data. □ 

The rest of the paper is devoted to the proof of Theorems [5] and [6] 



5. Formulation in the fixed domain 



In this section we transform the problem {P{n^)) to the fixed domain 51. 
Let us introduce the following notation: 

v,(t,x) := N^(x)vai,y(x)), voe(x) := NT(x)vo(y(x)), (t,x) e Q. 

Note that the definition of implies that divv^ = 0. The new function 



V, e L-(0,T;Wi;:;iJfi)) satisfies flN-iN-^atv, e [^L-{0,T;W'^'l,jn)) 
Ve(0, •) — voe and the equality 



f • 



foran0eL-(O,r;Wj':iiJfi)), (20) 



where the term on the right hand side is defined for (p e L^{0, T; Wq'^j^ (51)) 
by 



(g) N"^Ve : V(N"^0) - Ve ® Ve : V</> 
+ (C - 0N^^CN-^)ve • + isN-^f o y - f ) • 0I . (21) 



Here 



The left hand side of ()20p contains the perturbed term 0N~^N~^9(V£ due to 
the lack of a uniform estimate for Stv^. Similarly, the perturbed elliptic term 
0S(DeV£) : is present because of insufficient Lipschitz estimates, see Sec- 
tion [7] for more details. 
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Applying change of coordinates we get: 



I 



gN-^N-^ (ve(r) - voe) + [ \g (N^^CN^^v^ - N^^f o y) • | 



+ (NTS(©eVe) - Ve ® (N 




)) : V(N-^^) . (22) 



Now after all quantities and equations have been transformed to the fixed 
domain fi, we can analyze the limit e — > 0. 

6. Shape stability of weak solutions 

In this section wc prove that converges to v in certain sense. The result 
will be applied in the forthcoming sections. 

6.1. Uniform estimates 

Since the strong monotonicity of § and the Korn inequality hold uniformly for 
e 0, the energy inequality (HH) implies that ||ve||ioo(o^7^.L2(n,)), l|ve||^.(g_-r.wi.- 
and ||ve(r)||2.si^ is bounded uniformly with respect to e. Using the expansions 
we realize that 



where o(l)— ?>Oas£— >0. The same holds for the other norms, hence we derive 
the uniform estimates of v^: 



{vJ,>o is bounded in i°°(0, T; L2(17)) and in V'{0,T;Wlll.^{n)), (23) 

{ve(T)}e>o is bounded in L^i^). 
The Lebesgue-Sobolev interpolation inequality 



Using this information and (|6]) we can estimate all but the first terms in (|20p 
and thus 



I|ve(r)|i2,^., = llv,(r)ji2,o + o(i), 



||v,|i^,<C||Vv,||,||v,||2 



(24) 



together with yields that 



{ve}e>o is bounded in L^^{Q). 



(25) 



{f|N ^9tVe}£>o is bounded in 



(l'-(0,T;W^-,(17))) 
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6.2. Convergence 

The uniform bounds, the interpolation inequahty and the Aubin-Lions ar- 
gument give rise to the foUowing convergence: 

Ve-^v weakly-* in i°°(0,r;L2(f7)), 

weakly in L'-(0, T; W^;^; Jf})), 
weakly in L2'-(g), (26) 

Ve ^ V strongly in L^(Q), 1 < z < 2r, (27) 

V, (T) ^ v(T) weakly in (^2) , 

QE-'E-^^t^^e ^ 9tv weakly in (l'-(0, T; W^;^;^ (17)) 



N^§(I])eVe) ^ S(]D>v) weakly in L'' (Q, M^xa)^ 

- weakly in (^''(0, T; W^;:,.^ (17)) 



where v, S(Dv) satisfy the integral identity 

f (9tv, </))wi-'- (oi + / [s(]D)v) : - V ® V : V0 + Cv • J = / f • 

with the test functions 4> •= L'^{^iT;W^'^^^^{Q,)). Clearly v satisfies the initial 
condition v(0, •) = vq. To identify the weak limit S(Dv) we show that 

Dve ^- Dv strongly in L'^(g). (28) 

From jT]) we get: 



C5|lB(v,-v)|r;,Q< / (§(Dv,) 


- §(]D)v)) ; ©(ve - v) 


JQ 




= / (§(Ov,) 


- §(ID)eV.)) : ©(v. - v) 


Jq 




+ [ S(©eVe) 

Jq 


: (©(V, - V) - (NV(N-^(V, - ^)))sy,n) 


+ [ S(D,v,) 


: (NV(N-^(v,-v))),,y„ 


Jq 




- / §(©v) : 1 


])(v, - v) 


Jq 




= h+l2 + h 


+ /4. 



We immediately see that /4 — )- as e — >■ 0. Further ©eVe = Dv^ -I- eA^, where 
||A^||r,Q < C* uniformly w.r.t. e > 0, as follows from the uniform estimates for 
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Vg. Consequently it holds: 

Il= f [ -^S(PeVe + s(DVe - BeVe))ds : ©(v^ - v) 
/o CIS 



iQ 



1 

^'(BeVe + S(]D)V. - D£V£))<is(©Ve - DeVe) ; ©(Ve - v) 







<!1 ( §'(D,v, + ,s(Dv,-D,v,))ds|l_x_,Q|l£A;il.,Q||D(v,-v)|l,,Q 

< Ce, 







where we have also used (jSap . In the same spirit we obtain: 

h < Ce. 

The last term I3 can be expressed using ([20]) as 



/3 = - ^ (0N-iN-T9,v„ V, - v)^i,.^^ + 0(1), 

where o(l) — >■ as e 0, as follows from the available convergence and energy 
estimates. Integrating by parts with respect to time we get: 



h = II V0N-^(voe - vo)||^ - !| V0N-^(v,(r) - v(T))| 



2 



(flN-^N-^Stv, V, - v)^i.. + 0(1) = /31 + I32 + I33 + 0(1), 







where /31 — > 0, I32 < and 133 vanishes by the weak convergence of in 

L'-(0, T; Wl'^^in)). We have proved 

By means of the Vitali theorem we conclude that §(Dv) = §(]D)v) and con- 
sequently V = V is the weak solution to (P(f2)). As this solution is unique, the 
whole sequence {v^} converges to v. 

7. Existence of material derivative 

Our next task is to identify v as the limit of the sequence {u^}, where 

Ve V 
Ue := . 

e 

First we write down the system for the differences . Subtracting ([20)1 and ([5]) 
we find that Ue G L''-{Q,T;Wl['^.^{n)) satisfies sN-^N^^StUe G (l'-(0, T; WJl^j;^ 
Ue(0, •) = and the equality 



JQ 



e 

■ T 



1 1 f 

-Cue-0-(ve0Ue + Ue0v) : V0 = -/ (Ae, 0)^1.2 m) (29) 
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for all e L'^(0, T; Wl'^^^.^ {n)). The term A, G (i^^g, T; Wj;^;^ (17))) * on the 
right hand side is defined as follows: 



A, := AJ+A2+A3, 



A^ is given by (f2T|) . 
r{Al(l>)^.,. := /" (I - 0N-iN-^)atv . c^, 

^(A3, 0)^1,2^^ [ [n^§(©,v) : V(N-^0) - §(Dv) : 





Next we want to derive uniform bounds for {ug}. Therefore we refine the 
estimates of the previous section. In contrast to Section[Sl these estimates do not 
follow the structure of the problem, namely the elliptic term i 0(S(DeVe) — 

S(ID)ev)) : De(-) is bounded in (0, T; W^'f^ (rj)))* and the right hand side 

iA, in(i2(0,T;W^;2;,(r!)))*only. 

7.1. Lipschitz estimates for 

In this subsection we will investigate the Lipschitz continuity of the map 
e I— > A^. In particular, we arc going to show that 

e 

where Ag is defined in (|17ap . 

Indeed, due to ([T5|) and Corollary |4] it holds: 

^A^ -> (N' + N'^ - ItrN')5tV strongly in L^{Q). 

In order to pass to the limit in the terms -jAj and ^A^ we rewrite them in a 
convenient form: 

1 r"^ f /N^^ — I \ f — I 

'K,^)w'-l (U) - / ve®v, : V + / ve® V, : V 



A^ weakly in (^1^(0, T; Wl%^ (17))) * , (30) 



e J .Iq e 



£ 

6 



Jq ^ JQ £ ~—[ 



i j^{Al, 0)^1..^^ = ^ §(©v) : V (^^^J^0^ ^_i§(Pv) : V (N-^0) 



9 



3=^ 
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Employing we obtain for e 0: 



h ^ N'div (v (g) v) • 0, 

Jq 



div (v ® N'^v) • 4). 

Jq 

Similarly, from (P7)) it follows that 

i (N' -ItrN')Cv • 0, 



(ItrN'-N')f -c/). 



Analogously, 



h^- f S(Dv) : V (N'^0) 
Jq 



h I N'^§(Dv) : Oct). 



Rewriting /g carefully yields: 

/e = - / [(f • 0) o y + f o y . (0 - o y) - f . 0] 

—li-d,-- f {■4>+ f foy. 

^ ^Jq\Q'^ Jq 



QnQe 



0- </)oy 



/ [(f •0)divT + V0:f®T] =: / ((Vf )T, 0)^1,: 



where we have used the Hardy inequality and the following argument: 

|f.0| 



1 








[ f -0 




e 


JQ\Q, 


JQ\Q. 



distao 



<C!lf!l2,Q\QJ|V'All2,Q^0. 



Finally, taking into account the continuity of §' and Dv (as follows from 
Corollary 2]) , we have: 



/g = N^Miv (n'^ §'(]D)eV + .s(ID)v - Dev))ds — 



div (§'(Dv) (trN'Pv - (N'Vv)sa,„ +]D)(N'^v))) 



strongly in (^i'^(0, T; Wq'^j^ (fi)) j . Summing all terms up we arrive at ([30] 
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1.2. Estimates o/u^ 

Testing (p9)) by l[o,t)Ue and integrating by parts the time derivative and the 
convective terms we obtain: 



^ Jo Jn £ Jo Jn 

Next we apply (O, ([T5|). ([M)) . Holder's and Young's inequality to show that 
\\u,{t)\\l+ 1^ ||©u,||2<c(^ ||Vv,||2||Vu,||2||u,||2 + ||vo||^ 



for a. a. i e (0, T) and for t = T. Gronwall's and Korn's inequality and (|30 
then directly implies that 

{uJe>o is bounded in L°° {0,T;L^{n)) and in ^^(o, T; W^'^^^))^ 
{ue(r)}e>o is bounded in L'^{n). 
Since Dvg is bounded only in L''{Q), we have for e L*^ (0, T; W^'*^ (17)): 



3(S(eeVe)-§(ID)eV)) : = / g / S'(ID)eV+sD^(Ve-v))(isejUe : 

Jq Jo 

< C(||©v,||,,, + |l©v|l,)'-2||Bu,||2|l©0|| 2^ 



and thus 

■S(ID)^v^) -§(ID)ev) 



Consequently, due to (j^ . 



is bounded in L^--^ (Q). (31) 

:>0 



{oN-^N-^atUel^^^Q is bounded in (l^(0, T;Wo;£7 (17))) . (32) 



7. 3. Convergence of Ug 

We are ready to prove the existence of material derivatives. The estimates 
derived in the previous section imply that 



u. weakly-* in L°°(0,T;L'^(f2)), (33a) 

^OAiv ' 



weakly in i2(o,r;Wj:2.Jl])), (33b) 



weakly in L^((5), (33c) 
u.(T)^u(T) weakly in L2(f2), (33d) 

u. -> u strongly in L'^(Q), 1 < z < 4. (33e) 
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Since § G C^, the Nemytskii mapping associated to S' is continuous from L'^ to 
, hence from (1281) we deduce that 



-[ (0(S(ID)eVe)-§(ID'ev)),]D)£0) ^ /" (S'(Dv)Du, Be/.) (33f) 
^ Jo Jo 

for all (f> e i#^(0,r; W^'3^(J1)). Consequently, ([32]) yields: 

gN-'N-^dtu, - dtu (33g) 

weakly in ^L*^ (0, T; Wg'^j^'' (fi))^ . The initial conditions satisfy: 

Ue(0) ^ N'^vo + (Vvo)T = u(0) weakly in L^ifl), 
hence u = v is the unique weak solution to pT|. Theorem [5] has been proved. 

8. Shape gradient of J 

Let us decompose the fraction 

e 

into the sum Jf + J| where, in view of (|13p and 



J[:= / 0N-iN-^(u,(T)-u,(O))-| + / 



-rSflUeV.) -§(ID)v) , , , , 



and 



Jo £ Jq 



gN-^C v+i^ Cv 



-fe o y • ^ 



£ £ 
£ £ / 

/N"^ - 1 
+ (S(Dv) - V ® v) : V ( ^ 1 

From ^ and (|33l) it follows that 

Jl ^ Jv(v). 

Using similar arguments as in Section [73] we can show that 

J| ^ Je(T). 

The continuity of the map T i-> dJ{il; T) follows from the estimate (|18p . The- 
orem [S] is proved. 
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9. Conclusion 



The paper was devoted to the proof of shape differentiabihty for the prob- 
lem of time-dependent planar flow of an incompressible fluid with shear rate 
dependent viscosity. We have considered the so-called subcritical case when the 
weak solution is known to be unique and have bounded gradient. These strong 
properties were crucial in the proof of the existence of the material derivative, 
in particular for the well-posedness of the linearized system satisfied by v. For 
this reason it seems impossible to extend the results to 3 spatial dimensions and 
to r < 2, because no available theory can guarantee bounded gradient of the 
solution in such cases. 

The considered model itself has a limited practical applicability, mainly due 
to the homogeneous boundary condition. We believe that a generalization to 
other types of boundary conditions and shape functionals is possible to some 
extent, however to keep ideas clear the slightly artificial setting was chosen. 
Extension to the problems with non-trivial boundary conditions as well as ap- 
plication to numerical methods will be considered in the further research. 
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